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Steam pipelines applied in power units operate at high pressures and temperatures. In addition, to stress from the 
pipeline pressure also arise high thermal stresses in transient states such as start-up, shutdown or a load change of 
the power unit. Time-varying stresses are often the cause of the occurrence of fatigue cracks since the plastic de- 
formations appear at the stress concentration regions. To determine the transient temperature of the steam along 
the steam flow path and axisymmetric temperature distribution in the pipeline wall, a numerical model of pipeline 
heating was proposed. To determine the transient temperature of the steam and pipeline wall the finite volume 
method (FVM) was used Writing the energy conservation equations for control areas around all the nodes gives a 
system of ordinary differential equations with respect to time. The system of ordinary differential equations of the 
first order was solved by the Runge-Kutta method of the fourth order to give the time-temperature changes at the 
nodes lying in the area of the wall and steam. The steam pressure distribution along pipeline was determined from 
the solution of the momentum conservation equation. Based on the calculated temperature distribution, thermal 
stresses were determined. The friction factor was calculated using the correlations of Churchill and Haaland, 
which were proposed for pipes with a rough inner surface. To assess the accuracy of the proposed model, numer- 
ical calculations were also performed for the thin-walled pipe, and the results were compared to the exact analyt- 
ical solution. Comparison of the results shows that the accuracy of the proposed model of pipeline heating is very 
satisfactory. The paper presents examples of the determination of the transient temperature of the steam and the 


wall. 
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Introduction 


The high-temperature steam pipes applied in steam 
boilers are used to transport superheated live steam from 
the boiler to the turbine. 

During startups and shutdowns of the power units as 
well as during load changes high thermal stresses can 
occur in the pipeline wall and pipeline fittings. Particu- 
larly high stresses occur in T-pipes and Y-pipes. Also, 
thick-walled parts with complex shapes such as gate 


Nomenclature 


Dawid Taler: Dr. Science, Associate Professor 


valves and valve housings are exposed to large thermal 
loads. High and time variable thermal stresses may lead 
to premature damage to the pipeline in the form of cracks. 
Knowledge of the range of the stresses in critical com- 
ponents allows conducting a startup of the boiler in a way 
that provides a safe and long life of the boiler and turbine. 
Issues relating to the direct and inverse calculation and 
monitoring of thermal stresses in cylindrical components 
are the subjects of current research [1-2]. In this paper, 
transient temperature and thermal stress distributions in a 


a thermal diffusivity a=W/(c p), m’/s 


Ta nominal (design) temperature, °C or K 

t time, s 

time after which the steam temperature 
reaches the nominal temperature, s 


cn 


tpr transit time of the fluid particle, s 
XYZ Cartesian coordinates, m 
w fluid velocity, m/s 
Greek symbols 
a heat transfer coefficient, W/(m*-K) 
Ar radial step ,m 
At time step size, s 
A thermal conductivity, W/(m-K) 
E relative roughness, m 
dynamic viscosity, kg/(m' s) 
v kinematic viscosity, m/s 
Vr rate of change of temperature, K/s 
p density, kg/m? 
T time constant, s 
Ti time constant for wall, s 
Te time constant for steam, s 
Subscripts 
C steam 
in inner surface 
out outer surface 


mean 


A cross-section area, m? 
: specific heat capacity at constant presure, 
r J/(kg:K) 
specific heat capacity of the heat tube 
x material, J/(kg:K) 
din inner diameter, m 
l (x) modified Bessel functions of order n, 
L length of the pipeline, m 
m mass, kg 
m fluid mass flow rate, kg/s 
(n+1) number of nodes in the radial direction, 
n number of nodes in finite difference grid, 
NTU number of transfer units, 
Nu Nusselt number, 
p absolute pressure, Pa 
Pr Prandtl number, 
r radius,m 
Fin inner radius,m 
Foi outer radius,m 
Re Reynolds number, 
S wall thickness, m 
T temperature, °C or K 
fluid and tube wall temperature for linear 
T, Twr increase of the fluid temperature at the 
pipeline inlet, °C or K 
fluid and tube wall temperature for linear de- 
TmT wir crease of the fluid temperature at the pipeline 
inlet, °C or K 
fluid and tube wall temperature for a step 
Tm Twm change in fluid temperature at the pipeline 
inlet, °C or K 
T,, (0) mean temperature of the wall, °C or K 


pipeline connecting the boiler and turbine will be ana- 
lyzed. At first, the temperature of the fluid and pipeline 
wall temperature will be determined analytically using 
the superposition method. In the analytical solution, the 
physical properties of the fluid and metal are assumed to 
be independent of temperature. Also, a temperature 
difference across the thickness of the pipeline wall is 
neglected, i.e. the wall is modeled as an element with 
lumped heat capacity. Next, the finite volume method 
(FVM) will be used to calculate the transient temperature 
of the steam and pipeline wall. Thermal stresses will also 
be determined. Time and space temperature changes of 
the fluid and metal calculated using analytical and nu- 
merical methods will be compared. The pipeline heating 
with superheated steam will also be presented. 


w wall surface 


Simplified Analytical Model 


In a simplified mathematical model of heating or coo- 
ling of the pipeline, the following simplifying assump- 
tions were made: 

e the temperature of the fluid in the pipeline cross 
sections is constant and depends only on the time 
and axial coordinate, 

e the temperature drop across the wall thickness is 
negligible, 

e axial heat conduction in the pipeline wall has been 
omitted, 

e the outer surface of the pipeline is perfectly 
thermally insulated, 


e physical properties of the fluid and the wall material 
do not depend on the temperature and location. 


Equations for the steam and the pipeline wall are 
as follow: 
T 1 oT 
1+ =-(r-T,) (1) 
ðt NTU dz 
T 
Tw 7 = T T T, (2) 
ot 
Boundary conditions: 
T\ eT, + Ar +v OStSt, (3) 
T lal, t2t., 


where: 
T =1,+ AT + Vv, t, 
Initial condition 
T |= Tp 
T, liso = T 
The superposition method was applied to solve Eqs. 
(1-4). Formulas for the transient temperature of the liquid 
and the wall for the abrupt and a linear increase in fluid 
temperature at the inlet of the pipeline are given in [3]. 
The solution of the boundary-initial value problem 
given by Eqs.(1-4) is: 
The temperature of the steam 
T=T, (+p), Ostet, > ta Zt, ©) 


T=T,(t)+ Ty (tbe) + TiO), t Sten (6) 
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Fig. 1 Time changes of the fluid temperature at the pipeline 
inlet- boundary conditions for solution components I, 
II, and II in the superposition method 
where: 
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The temperature of the pipeline wall 
Ty pon Tom), OSCE Tila ty (10) 


Tg Qed, w(t- fen) + Tun (11) 
where: 
Tyr =T, -vpr Ë ”) KE 1) + Io (2Vén) | (12) 
go Jt S 
Ton (ten = Ty + vere ee [U (ém) 
+o (2/ém ) | (13) 
ee on Z tpr 
Ton = 19+ AT e ve) lu (67 )+1,(2 En) | (14) 
t2 tpp lon È bpr 


The following designations in formulas 1-14 were 
adopted [3]: 


-eln _ $2}, (2/é7) 


v ; U (é, ) 
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pa ZNTU t-t 
L’ 7 ‘ > 
t=t,, =,» 
m — tp =Z NTU Tos 
w 
NTU = As ze aA, : 
MC ye PWA, Coc 
E Cpe m E p-LA, Cpe 


P 

‘ ad, aA, 
The analytical solution given by Eqs. (5-14) will be 

used for comparison with the numerical solution pro- 


posed in this paper. 
Numerical Model 


A scheme of the pipeline linking the boiler with a 
turbine in power unit with capacity of 120MW is 
depicted in Fig. 2. 

The governing equations for the steam are 
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mass conservation equations 


ep __1 h (15) 
ôt AO 
momentum conservation equations 
s aD) 
an =— o2 -a{ 2+ pgsino 
Ot 0z\ pA Oz 
Ea (16) 
„Z me 
d 2 p £ 
energy conservation equations 
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ot pAda ôt pAðz 
(17) 
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Fig. 2 Pipeline connecting the boiler and turbine 


Transient heat conduction equation for the pipeline 
wall is as follows 

ôT, = V- A, (7) VT, 18 

PS pw at = [ al w) p] ( ) 

Taking into account small changes in pressure along 

the length of the pipeline, a small steam density change 


can be assumed that the derivative Op/Ot is equals to 
zero. The assumption 09/0t=0 in Eq. (15) implies 
that mass flow rate m is constant over the tube length, 
Le. 
m = const (19) 
The velocity of the steam in each cross section of the 
pipeline can be calculated from the formula 
we (20) 
1 Pdin 


Considering that Om / Ot = 0 , the momentum con- 
servation equation (17) reduces to 


dp ing Že" 
dz p dz PRR d. 2 


in 


(21) 


Using Eq. (21) the steam pressure can be determined 
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along the flow path. 

Neglecting the thermal expansion of the steam (B = 0), 
the heat generation due to friction and the axial thermal 
conduction in the steam, the energy conservation equa- 
tion (17) for the steam is simplified to form 

T n OT a(T—-T,,)U,; 
Pe, ôT m ôT __2(T-T,)Un (22) 
ôt pA@z A 

The heat conduction equation (18) in cylindrical coor- 

dinate system is 


ôT, -1%] (r | 


Pw pw a, = r 
ot ror ; or (23) 
(a) E 
+— T, )—— 
Ha w) a) 


System of partial differential equations (22) - (23) is 
subject to the following initial and boundary conditions 


T |o=T, (24) 

Ti l=0 = Tio (25) 

T l= f(O) (26) 

E |. =a(F=7,) (27) 

ðr 
ay |, = 9 (28) 
Or u 

T 

n ijt 29) 
z 

1 aL (30) 
Oz 


In addition to the conditions (24) - (30), the pressure 
and mass flow rate of the steam are known at the inlet of 
the pipeline. 

The initial -boundary value problem (24)-(30) was 
solved using the finite volume method (FVM)[4-6]. 

First computational domain, i.e. the wall and the area 
occupied by the steam was divided into finite volumes. 


(Fig.3.) 
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Fig. 3 Division of the computational domain into finite vol- 
umes. 


For each control volume lying both in the wall and the 
steam, energy balance equations were formed. For exa- 
mple, the energy balance equation for node i is set for 


control volume located in the computational area of the 
wall (Fig.4.) 
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Fig. 4 Wall temperature at the node i and adjacent nodes i-1, 
i+1, i-n-1, itn+1. 


The energy balance equation for the node i is 
aT,, 
Pw (7, Wi Je Cow (T, Jafra- r? )Az P 


+27r AZ Ay (Tn) (Pan) i fe, (1) 
ta(r2-72) Ay (Ty,)+ = ae Mea Ta 
ale ea Tn To 


Transformation of Eq.(31) gives 


dT, a(Ty,) | r; 2y (Zn J+ Aw (Tu, ) a, 
dt Ay (Ty, } Ar Pin mW 


a pag a a 
Avy a )+ Avy ae ) (7, T ) 

2 (Az) Witn+1 Wi 
Any (G j Ay (Bis ) (T, T ) 

2 (Az) Wi-n-1 Wi 


In a similar manner, the heat balance equation for the 
i-th control volume located in the region occupied by 
steam can be set 


aT, = m T4 —T; a(T; U, Tau-i 
dt p(T;)A Az = Ap(T;)c, (7; 
m) Mell 2. 
miaii wa 
i(n-1)+1 (i-1)(n+1)+1 i= 1,...,m 
2 | 


The heat transfer coefficient a on inner surface of the 
pipeline was determined using Gnielinski’s correlation [7] 


Š Re-1000) Pr 2/3 

Nu= a - (4) | (34) 
1+129( £) (Pr°°-1) 
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where the friction factor € is given by the Colebrook 
correlation 


é =[1.8log(Re)—1.15] * (35) 

The Reynolds number, Prandtl number, and Nusselt 
number are defined as follows 

Pr = <2", Nu= oes 
u A A 

After the formulation of all balance equations for the 
wall and steam, a system of ordinary differential equa- 
tions for node temperatures is obtained. The number of 
equations for the nodes lying in the wall area is 
(m+1)(n+1) while for the steam is (m+1). 

The system of ordinary differential equations was 
solved by the Runge-Kutta method of the fourth order. 

To ensure the stability of determining the wall and 
steam temperature, the Fourier stability condition for the 
wall and Courant-Friedrichs-Levy condition for the steam 
should be satisfied. The smallest allowable time step 
Atmax results from the Courant-Friedrichs-Levy condition 

w;At 
Az 

Knowing the transient temperature distribution in the 
wall thermal stresses can be determined. 

Considering that the axial component of the tempera- 
ture gradient OT,,/0z is very little only radial tempera- 
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Re= (36) 


<1, i=1,...,m+1 (37) 


ture drop in the wall of the pipeline is taken into account. 
Assuming that the ends of the pipeline are free to 
lengthen the thermal stress components are given by the 
following formulas [8] 
Epr 
0, = 1-2 || TŒ -TC,t 38 
F ate :|[To- Te) | 69 


EBr in r 
o= E rosis thien- zre] (39) 


EBr c= 
7, =| TY) -2TG,t 40 
o: =a p [TO-2769] (40) 
where the mean wall temperature T(t) 


2 


where: 
r=1+Ar/2 
r=1r,+(-2)Ar; 
Ty =r a tAr/2 


n n+l 


i=3,... n+l1 


Radial stresses o, are equal to zero on the inner and 
outer pipeline surface. Circumferential o, and axial 
stresses ©, are equal on those surfaces. 


Results 


This paper presents the results of calculations of the 
pipeline connecting the boiler OP-380 with a steam tur- 
bine. 

The steam pipeline is made of low alloy steel 13HMF 
(C-0.18%, Mn-0.40%, Si-0.35%, Pmax-0.040%, Smax- 
0.040%, Cumax-0.25%, Cr-0.60%, Nimax-0.30%, Mo- 
0.65%, Almax-0.020%). The main dimensions are outer 
diameter d,,,;= 0.324 m, the wall thickness s = 0.04 m, 
and length L = 45 m. 

First, the results obtained by the numerical method 
proposed in this paper will be compared with the analyt- 
ical solution. The comparison will be carried out for con- 
stant physical properties of steam and steel because of 
constant physical properties assumed in analytical solu- 
tion. The following constant physical properties have 
been used in calculations: 

Aw = 44.68 W/(m-K), py = 7766 kg/m’, cy = 545.3 
J(kg-K), a = 0.0133 m’/s, Cpe = 2672.3 J/(kg-K), po= 
39.8 kg/m’, 4. = 0.0837 W/(m-K), v = 0.7755-10° m’/s, 

The initial temperature of the fluid and the pipeline 
wall was equal to Tọ = 20°C. Initial temperature jump 
was AT = 100K. The steam temperature increased at a 
constant rate vr= 10 K/min = 1/6 K/s in the range from 0 
to ten = 2519.9 s. Steam velocity in the pipeline is con- 
stant: w.= 56.57 m/s (Re = 1.77989410’). Heat transfer 
coefficient was calculated from the Dittus-Boelter for- 
mula: 


Nu = 0.023Re°* Pr’ (43) 

Given that the Reynolds number is Re = 1.779894x10’, 
and the Prandtl number Pr = 0.985, the heat transfer 
coefficient calculated using Eq.43 is: 

a = 4959.2 W/(m’-K). The number of transfer units is 
equal to NTU = 0.6 and the time constants for the fluid 
and pipeline wall are: = 1.31 s, t= 39.76 s, respec- 
tively. 

Comparing the temperature of the fluid and the pipe- 
line wall calculated using the analytical formulas and the 
proposed numerical method is depicted in Fig.5. From 
the analysis of the results shown in Fig. 5 it can be seen 
that in the initial stage of pipeline heating, the mean wall 
temperatures calculated from the analytical formula and 
by using the finite volume method differ significantly. In 


a further heating phase the pipeline, average temperature 
of the pipeline wall calculated from the analytical formu- 
la is similar to the internal surface temperature calculated 
using the numerical method developed. 


600 
(a) 
400 ie 
y ; 
E ff 
200 Se Gates T, z=0m 
5” T, z=45m (anal. sol.) 
Tw, z=45m (anal. sol.) 
©---O--OT, z=45m (n.model) 
* -X - X Twm, z=45m (n.model) 
——__— Tw, r=r,,, Z=45m (n.model) 
0 
0 1000 2000 3000 
ts 


T, z=45m (anal. sol.) 
Tw, z=45m (anal. sol.) 

| _|Q---O--OT, z=45m (n.model) 

X -- X Twm, z=45m (n.model) 
-—— Tw, f=1,,, Z=45m (n.model) 


Fig. 5 Comparing the temperature of the fluid and the pipe- 
line wall calculated using the analytical formulas and 
the proposed numerical method; a) the time interval 
from t = 0 to t = 3000s, b) the time interval from t = 0 
to t= 500s 


Differences between the temperatures of fluid calcu- 
lated using analytical and numerical methods are slightly 
different. 

The numerical method must be used to calculate the 
thermal stresses in the wall of the pipeline. 

Test calculations of the pipeline heating were per- 
formed taking into account the temperature-dependent 
physical properties of the fluid and the pipeline material. 

Thermo-physical properties of steel are a function of 
temperature. 


A = 48.495 +0.0012T —6.46x10°T7 


ro (44) 
+4.175x10°T 


a =1.341x10° —3.452x10°T 


—3.193x107!! T? +2.853x10°4 73 


where the thermal conductivity 2 is in W/(m-K), the 
thermal diffusivity a in m’/s, and temperature T in °C. 

Calculating the transient temperature of the fluid and 
pipeline wall was carried out for a different number of 
finite volumes across the thickness of pipe wall. 

Inspection of the results shown in Fig.6 illustrates that 
even at four finite volumes, i.e. for the five nodes evenly 
distributed over the wall thickness satisfactory accuracy 
of the calculations is obtained. 

Almost identical results are obtained for the division 
of the pipe wall into nine (10 nodes) or nineteen finite 
volumes (20 nodes) as for the division into four finite 
volumes (5 nodes) Fig.6. 

The temperature of the wall in the radial direction is 
computed in five nodes evenly spaced (n=4). The number 
of nodes in the axial direction is equal to m+1 = 21. 

The pressure is p = 13.9 MPa. The mass flow rate of 
the steam is equal tom =105.55 kg/s. The initial tem- 
perature of the pipeline and steam is Twy = 20°C. At time 
t > 0, the temperature of the steam pipeline inlet increas- 
es abruptly to a constant temperature at Ty>=540°C. 

Selected modelling results are shown in Figures 7-11. 
The steam temperature as a function of time at nodes No. 
2 (z = 2.25 m), No.10 (z = 20.25 m), and No. 20 (z = 
42.75 m) are displayed in Fig. 7. 
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Fig. 6 The temperature of the inner and outer surface of the 
pipeline at its inlet (z = 0) as a function of time ¢ for 
various numbers of control volumes n across the 
thickness of the pipeline wall. 


The analysis of the results shown in Fig. 7 shows that 
after about 600 s steam temperature reaches a steady 
state. 

Fig.8 illustrates steam temperature changes over the 
length of the pipeline at time points 10 s, 60 s, 240 s, and 
600 s. 


The pipeline wall temperature in two different cross- 
sections as a function of time depicts Fig.9. The wall 
temperature is shown in five uniformly spaced nodes. 
The nodes 11 and 101 are located on the inner surface 
while the nodes 15 and 105 lie of the outer pipeline sur- 
face. 
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Fig.7 Steam temperature as a function of time at nodes No. 2 
(z = 2.25 m), No.10 (z = 20.25 m), and No. 20 (z = 


42.75 m) 
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Fig.8 Steam temperature changes over the length of the pipe- 


line at time points 10 s , 60 s, 240 s, and 600s. 


At the beginning of the heating process, the tempera- 
ture differences over the thickness of the wall are large 
but rapidly decrease in time (Fig.9.). 

Circumferential stresses on the inner and outer surface 
of the pipeline as a function of time presents in Fig.10. 

The analysis of the results presented in Fig.10 shows 
that the inner surface is exposed to high compressive 
stresses while on the outer surface occur substantially 
lower tensile stresses. 

High thermal stresses in the pipeline are caused by a 
high jump of 520 K in steam temperature at time /=0. 

The distribution of the circumferential thermal stress 
over the wall thickness is shown in Fig.11. 
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Fig. 9 Wall temperature in two cross-sections as a function of 
time a) z = 2.25 m, b) z= 42.75 m. 
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Fig. 10 Circumferential thermal stress at the inner (nodes 11 and 
96) and outer (nodes 15 and 100) surface of the pipeline 
as a function of time; (a) z = 2.25 m, (b) z = 42.75 m. 
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Fig. 11 Circumferential thermal stress as a function of radial 
coordinate and time (a) z = 2.25 m, (b) z = 42.75 m. 


The stresses are compressive near the inner surface. In 
the vicinity of the outer surface, the stresses are tensile. 

The largest absolute value of stress occurs at the be- 
ginning of heating and tend to zero over time. 


Summary 


The mathematical model of the steam pipeline heating 
developed in the paper allows determining the tempera- 
ture of the steam and the pipeline wall as a function of 
position and time. A comparison of the transient steam 
and pipeline wall temperature obtained from the analyti- 
cal formulas with the results obtained using the proposed 
numerical model was conducted. This comparison indi- 
cates that the analytical formulas can only be used to ca- 
lculate the temperature of the steam. It is possible to 
determine transient thermal stresses caused by the tem- 
perature difference across the wall thickness, using the 
model developed in the paper. Examples of calculations 
of the steam temperature, wall temperature and circu- 
mferential thermal stresses on the inner and outer surface 
of the pipeline have been presented. The calculation tests 
performed in the paper shows that the developed mathe- 
matical model can be used to simulate the actual pipeline 
heating or cooling in a power plant. 
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